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The Initial Core Mass Function due to Ambipolar Diffusion in 
Molecular Clouds 
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ABSTRACT 

We show that the ambipolar-diffusion-initiated fragmentation of molecular clouds leads sim- 
ply and naturally to an initial core mass function (CMF) which is very similar to the initial 
stellar mass function (IMF) and in excellent agreement with existing observations. This agree- 
ment is robust provided that the three (input) free parameters remain within their range of 
values suggested by observations. Other, observationally testable, predictions are made. 

Key words: diffusion — ISM: clouds — ISM: individual (Orion) — magnetic fields — stars: 
formation — stars: mass function 
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1 INTRODUCTION 

Observations reveal that the initial stellar mass function (IMF) may 
be approximated by a power law over two decades in mass, with 
indications of a turnover below ~ 0.5 Mq and an upper mass limit 
~ 50 Mq. In this mass range, the number of stars dN per stellar 
mass interval dm* may be represented by d N/dm* oc m* a , with 
ol — 2.5 i 0.3. This range for a includes the ISalpeteiU 19551) IMF 
(a — 2.35) and also allows for a slig ht steepening at the h igh- 

see 
and 



mass end, as in th e Miller-Scalo IMF {Mil ler & Scalo 
reviews by| .Kroupa 2002, Bonne ll. Larson & Zinnecker 
iMcKee & Ostrikeill200l § 3.3). 
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Submillimeter observations of molecular-cloud cores show 
that, if the initial core mass function (CMF) is written as 
dN/dm c oc rxi^ a , then a ~ 2 — 2.5 for core masses m c > 
0.5 Mq, with a = 2.35 being typical, while a ~ 1.5 for 
m c ;$ 0.5 Mq (see review bv lWard-Thompson et alj|200l § 5).Q 
The fact that the CMF is very similar to the observed IMF suggests 
that the early fragment ation process may determine the stellar mas s 
spectrum as well (e.g.. lSimpson. Nutter & Ward-T hompson 2008). 

Since ambipolar diffusion is an unavoidable process in self- 
gravitating, weakly-ionized, magnetic systems, such as molec- 
ular clouds, and since the ambipolar-diffusion theory has had 
many of its predictions confirmed by observations (e.g., s ee re- 
views bv iMouschoviasI Il987l Il99fj : Iciolek & Basi] l200d : also, 
iMouschovias. Tassis & Kunal2006l § 4.2.1), one would expect that 
the IMF should also be a natural consequence of the ambipolar- 
diffusion-initiated, single-stage fragmentation of molecular clouds. 



1 Ijohnstone et all l2000l) caution that the flattening below 0.5 Mq in p- 
Ophiuchi m ay be from incompleteness due to limited sensitivity around 
~ 0.4 Mm: lMotte. Andre & Ner] ll998l) also do not find a mass turnover 
down to their completeness limit ~ 0.1 Mfl. However, SCUBA data from 
Orion show a mass turnover ~ 1 Mq I Nutter & Ward- Thompsonll2007l) . 
well above the completeness limit, suggesting that the turnover is not a se- 
lection effect. 



IMouschoviasI jl987t Il99lh showed analytically, based on a 
comparison of three natural lengthscales in molecular clouds, 
that ambipolar-diffusion-initiated fragmentation selects protostel- 
lar masses in th e appro xim ate range 1 — 30 Mq . More rec ently, 
iNakamura &~Lil J2008h and iBasu. Ciolek & WursteJ J2009h have 
extracted CMFs from numerical simulations of fragmentation in 
weakly-ionized, isothermal, magnetic molecular clouds that are in 
rough agreement with observations. 

In this paper we show that the CMF is a natural conse- 
quence of ambipolar-diffusion-initiated fragmentation in molecu- 
lar clouds. 



2 PHYSICAL CONSIDERATIONS 

Both theoretical calculations of self-gravitating, mag netic clouds 
embedded in a hot and tenuous ex ternal medium ( IMouschoviasI 
1 19761: [Fielder & Mouschovias! 1 19931) and analysis of observations 
dBasdl2000l : iTassi j|2007h~ have shown that such clouds are oblate 
and preferentially flattened along the magnetic field. We thus con- 
sider a weakly-ionized, isothermal, magnetic molecular cloud of 
temperature T, embedded in an external medium of constant pres- 
sure Pcxt, and threaded by a magnetic field B — B(x, y, t)z that 
has negligible variation along the 2— axis. We denote the neutral 
mass column density along the field lines by a n (x, y, t), the half- 
thickness by Z(x, y, t) [= a n /2p n , where p n (x, y, t) is the neutral 
mass density], and the ion number density by n\ (x, y, t) . Force bal- 
ance along the field lines is then expressed by 



fcxt + 2^ a ^ ' 



(1) 



where C = {k-^T /m^) 1 ^ 2 is the isothermal sound speed and m n 
is the mean mass of a neutral particle (= 2.33 amu for an H2 gas 
with a 20% He abundance by number). 

The effect of the external pressure on the cloud (of noncon- 
stant half-thickness Z) can be represented as an increase of its 
isothermal sound speed from C to C e ff: 
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C!s = n 2 [3Pext + Qr/2)gog] 
G 2 " 2 an [Pext + (7r/2)G(Tg] 2 



1 + 3P e: 



(2) 

(1+Pcxt) 2 

(Mortor3 [l99ll ; also lCiolek & Mouschoviasll 19931 eq. 28b), where 
the dimensionless quantity P cxt = 2P ext /irGo 2 is the ratio of 
the external pressure and the vertical self-gravitational stress. If 
Pcxt -C 1, the cloud is held together by its self-gravity alone and 
Ceff — * C 2 > whereas if P cxt 3> 1, the cloud is confined by a large 
external pressure, which can significantly affect the cloud's evolu- 
tion. Typical self-gravitating molecular clouds have P cx t — 0.1 
(see Mouschovias & Morton 1992, §2.1 for explanation). Using 
equation (O, e quation ([TJl may be solved fo r the column density 
and written as dCiolek & Mouschovia j[l993l eq. 72) 



1.15 x 10 

(1 + Pcxt) 



— ( n ° ) 
1/2 V10 4 cm- 3 / 



1/2 



10 K 



1/2 



gem 



(3) 



where n n {x, y, t) — p n /m n is the number density of neutrals. 
The (elastic) collision time of a neutral particle with ions is 



1.4 



ra\ ) ni{aw)m 2 



(4) 



The ion mass is m\ = 25 amu, the mass of the typi- 
cal atomic (Na + , Mg + ) and molecular (HCO + ) ion species in 
clouds; (<7Ui}iH 2 is the neutral-ion collision rate, equal to 1.69 x 
10 -9 cm 3 s _1 for HCO + — H2 collisions and almost identical to 
this value for Na + — H2 and Mg + — H2 collisions. The fa ctor 1.4 
accounts for the inertia of helium (see Mouschovias 1996). We as- 
sume the canonical relation between m and n n : 



K 



( n " ) 
V 10 5 cm" 3 / 



1/2 



with JC ~ 3 x 10~ 3 cm -3 . 



(5) 



The ionization-equilibrium parameter K, oc (f/aar) may take 
on a range of values, depending on the local ionizatio n rate C and 
the dissociative recombination rate a^r — 10 " 6 s" 1 dLanged 19851 ; 
Dalga rnol 1987b . The ionization rate in the dense interiors of molec- 
ular clouds is primarily due to cosmic rays, for which observations 
imply a typical value £ ~ 5 x 10~ 17 s -1 , with an approx imate 
factor of 2 uncertainty (e.g., Payn e. Salpeter & Terzianll 19841) . 

The thermal critical lengthscale is a natural lengthscale in 
molecular clouds and is given by 



Ga n 



(6a) 



0.15 



(1 + 3P cxt ) 
(1 + Pext) 3 / 2 



10 4 cm' 3 N 1/2 



T 
10 K 



1/2 



pc; 



(6b) 



it is similar to the Jeans lengthscale. lMouschoviasI dl987l[l99ll) ar- 
gued that the relative magnitudes of this and two additional length- 
scales (the Alfven lengthscale and the magnetic critical length- 
scale) play a crucial role in initiating fragmentation in molecular 
clouds and in determining the masses of the resulting fragments. 
However, no mass spectrum was calculated. The total mass con- 
tained within a radius equal to the initial At is 



rriT 



G 2 a n 



(7a) 



3.91 



(l + 3Pcxt) 2 /!0 4 cnT 



(1 + Pcxt) 7 / 2 



1/2 



10 K 



3/2 



Mr: 



(7b) 



We refer to this mass as the thermal critical mass. 



The mass-to-flux ratio of a cloud relative to its critical value 
for collapse, i.e., 



CTn/B 



l/(63G) 1 /2 l/(63G) 1 /2 ' (8) 

which is a result of nonlinear calculations dMouschovias & S pitzer 
1 19761) . is a measure of the importance of magnetic fields in the sup- 
port and evolution of molecular clouds. Linear stability analysis 
of a uniform, self-gravitating disk threaded by a uniform magnetic 
field perpendicular to th e plane of the disk gives a critical mass- 
to-flux ratio, 1/2ttG 1/2 dNakano & Nakamura|[l978l) . 25% greater 
than the nonlinear result. 



3 AMBIPOLAR-DIFFUSION-INITIATED 
FRAGMENTATION 

3.1 The Fragmentation Lengthscale 

The three lengthscales mentioned above can be recovered as limit- 
ing cases of a more general result. If we take the model cloud dis- 
cussed in § 2 to be a uniform, static, background state with small- 
amplitude perturbations parallel to the equatorial plane of the cloud 
(which introduce sinusoidal ripples on the cloud's surfaces), the 
dispersion relation governing suc h perturbations is given by (see 
Martodll99lt ICiolek & Basufl2006l) 



u) 2 + 27rGfc(j n ,o — k 2 C 2 tf 



uj + ik 2 v% fi T ni fi I 1 + 



kZo) 



- ^fc 2 fi.o ( 1 + 



kZ 



(9) 



where u) and k are the frequency and wavenumber of the perturba- 
tions, and va = B /(Anp-n) 1 / 2 is the Alfven speed. Quantities with 
the subscript "0" refer to the zeroth-order (background) equilibrium 
state. This dispersion relation contains three dimensionless free pa- 
rameters: the initial mass-to-flux ratio /xo in units of the critical 
value for collapse of a uniform thin disk, 1/27TG 1 / 2 ; the ratio P ex t 
of the external pressure and the initial vertical self-gravitational 
stress, (ty/2)Go 2 ; and the (dimensionless) initial neutral-ion col- 
lision time, 



Tni.O 



T"ni,0 



ix, 

0.241 



(1 + P c: 



U/2 



3 x 10 -3 cm -3 
K. 



(10a) 



(10b) 



For the canonical value P cx t = 0.1, the unit of time (G/27rGa n ,o) 
is equal to 0.385rff,o, where ts,o = (37r/32Gp n ,o) 1,/2 is the 
(spherical) free-fall time at the density p n ,o of the background state. 
If we take the cosmic-ray ionization rate to be in the observation- 
ally reasonable range £ = 2.5 x 10 — 17 — 10~ 16 s _1 (e.g., Dalgarno 
l2006h . it follows that f n i,o = 0.162 - 0.325. 

A mode is gravitationally unstable if the imaginary part of its 
frequency u is positive. We denote the wavelength (= 2n/k) of 
maximum growth rate obtained from the dispersion relation ® by 
Ag^ax, and we refer to half of this wavelength as the fragmenta- 
tion lengthscale, Af r = A g . max /2. We identify this with the initial 
radius of a molecular-cloud fragment. 

Figure [T] exhibits the ratio Af r /AT as a function of /iq for the 



2 One can verify this by choosing natural units as follows: [velocity] = C; 
[length] = C 2 /2 7rG<7 nj o; [column density] = cr n! o; and [magnetic field] = 
Bo; hence, [time] = C/27rG<7 n ,o; and [mass density] = 2nGa 2 l q/C 2 . 
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Figure 1. Ratio of the fragmentation lengthscale Af r and the thermal critical 
lengthscale A-p as a function of the initial mass-to-flux ratio fio (normalized 
to the critical value for collapse) for P ox t = 0.1. The different curves 
correspond to T n io = 0.162 (dashed line), 0.230 (solid line), and 0.325 
(dotted line). 



Figure 2. Core mass m c (normalized to the initial thermal critical mass 
m.T o) as a function of the initial mass-to-flux ratio fio (normalized to the 
critical value for collapse) for P ox t = 0.1. The different curves correspond 
to ? nii0 = 0.162 (dashed line), 0.230 (solid line), and 0.325 (dotted line). 



fiducial case P cxt = 0.1 and T n i,o = 0.230 (solid line). For very 
subcritical or very supercritical fragments, Af r — > At. In the former 
case, ambipolar diffusion allows the neutrals to contract quasistat- 
ically under their own self-gravity through the field lines until a 
thermally and magnetically unstable core forms. In the latter case, 
the magnetic field is dominated by the self-gravity of the fragment 
from the outset and a supercritical core contracts dynamically on 
a magnetically-diluted free-fall timescale. For marginally critical 
fragments, instability occurs via a hybrid mode. An important con- 
sequence is a dramatic increase in the fragmentation lengthscale, 
depending principally on the value of T n i,o, which is in turn deter- 

1/2 

mined by the ionization-equilibrium parameter K, oc rii/n n (see 
eqs. 11 Obi and (5))- We therefore also plot the fragmentation length- 
scale for r n i,o = 0.162 (dashed line) and 0.325 (dotted line). The 
greater (smaller) the ionization-equilibrium parameter, the greater 
(smaller) the maximum dimensionless fragmentation lengthscale. 
Since we are not concerned here with pressure-confined molec- 
ular clouds, we do not vary the paramete r P ex t ; the e ffect of a 
large P ex t on the sol ution is discussed bv lMortor] i 1 99 lh (see also 
Ciole k & Bas32006h . 



3.2 The Core Mass 

We define "cores" as self-gravitating, molecular-cloud fragments 
that have become magnetically supercritical by some means (e.g., 
ambipolar diffusion). This theoretical definition is complementary 
to the observational definition, which refers to a "subset of starless 
cores which are gravitationally bound and hence are expected to 
partic ipate in the star formation process" JWard -Thompson et al.l 
5). The core mass is therefore the total mass contained 
within a radius equal to the fragmentation lengthscale Af r , cr , 



where the subscript "cr" denotes quantities evaluated at the mag- 
netically critical state for collapse. 

For initially subcritical fragments (i.e., fio < 1), B 
depends weakly on density during the fragmentation stage 
dFielder & Mouschoviasl 1 1993L Fig. 9c). The following relations 
then hold: 

O"n,0 _ Bq B, 
TkT ~ 27TGV2 



(12) 



\2 



(ID 



2ttG 1 /2 - 

(see also iMouschoviasll 1 99ot eq. 5). For fragments that begin their 
lives as critical (i.e., /io = 1), &n,o = o"n,cr-U 

We rewrite the core mass i ll It by using equation d 12b to elim- 
inate (J n , cr in favor of the initial (dimensionless) mass-to-flux ratio 
fio and the initial column density cr n ,o, and then we use equations 
d6ab and fTat to introduce the initial thermal critical mass tot,o, to 
find that 

m c ~ min(/x , 1) (J^j m T,o = f(Ho, T n ;,o) Wt.o • (13) 

(The ratio A&-/At is independent of the stage of evolution of the 
core; hence, the subscript "cr" has been removed from these length- 
scales in eq. 1131 .) Any deviation in m c from mT,o is due solely to 
the presence of the magnetic field, whose effect is communicated 
mainly through the initial mass-to-flux ratio and, near the critical 
value of the mass-to-flux ratio, by the neutral-ion collision time as 
well (see Fig. 2). 

The dependence of the core mass m c (normalized to the ini- 
tial thermal critical mass tot,o) on fio is shown in Figure [2] for 
P cxt = 0.1 and Tni.0 = 0.162 (dashed line), 0.230 (solid line), 



J The column density enhancement can be obtained rigorously by calcu- 
lating the eigenvectors of equation (5). We have followed this approach to 
obtain the relations plotted in Figures[2]-H but we present only the approx- 
imate equation 1 1 2t in the interest of brevity and simplicity; the two results 
differ by at most 6%. 
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Figure 3. Predicted CMF (normalized to the initial thermal critical mass 
m T o) f° r Pcxt =0.1 and r n i o = 0.162 (dashed line), 0.230 (solid line), 
and 0.325 (dotted line). T he corresponding power-law expon ent a is also 
shown. The data are from iNutter & W ard- Thompson (2007) and refer to 
starless cores in Orion; see text. 



and 0.325 (dotted line). For Ho < 0.5, m c ~ "1t,o> whereas 
for no 3> 1, m c ~ rriT.o- In-between these two extremes, the core 
mass increases significantly due to the action of gravitationally- 
driven ambipolar diffusion. 



3.3 The Initial Core Mass Function 

In the absence of magnetic fields, the core mass is equal to the 
thermal critical mass given by equation {Taj, i.e., m c .B=o = 
rriT.o oc n~j' 2 T 3 ' 2 . Given that the deep interiors of molecu- 
lar clouds (where stars form) are observed to be approximately 
isothermal, the width of the CMF is then determined by the width 
of the initial probability density function (PDF) of number den- 
sities. In order to generate cores spanning almost two decades in 
mass, one therefore requires an initial density PDF spanning al- 
most four decades. By contrast, in the presence of magnetic fields, 
m c oc f{fio, Tni,o)mc,B=o- Therefore, one may obtain a spectrum 
of masses due to changes in the local mass-to-flux ratio. In fact, 
one only requires a one-decade spread in initial mass-to-flux ratios 
(0.1 < fio < 1) to generate a two-decade spread in core masses 
(see Fig.O. 

To derive the CMF predicted by the theory of ambipolar- 
diffusion-initiated star formation, we make the following assump- 
tions: (i) All protostellar fragments form by ambipolar diffusion 
with initially subcritical or critical mass-to-flux ratios (i.e., ho < 
1). (ii) The PDF of initial mass-to-fiux ratios is broad, in the sense 
that no ~ 1 is almost as likely as ho — 0.6 (we explain below the 
significance of the value Ho ~ 0.6). (iii) The linear analysis we use 
to calculate core sizes and masses is adequate to predict the frag- 
mentation properties of clouds even in a fully nonlinear stage of 
evolution. This is supported by comparisons of numerical simula- 
tions and predictions from the linear theory dCiolek & B asu 2006; 
iBasu, Ciolek & Wurstej|2009h . 

We then calculate the PDF of core masses V(fh c ) generated 
from an assumed PDF of initial mass-to-flux ratios V(ho) (for a 
fixed r n i,o): 



V(m c )dlnm c = Ho~P(fJ,o 



dim 



d In m c 



(14) 



where ra c 



d In no 

m c /m,T,Q- Making the reasonable assumption that 



Figure 4. High-mass end slope a versus the maximum core mass m c , max 
(normalized to the initial thermal critical mass mx,o) for P ox t = 0.1. Each 
point represents a different value of r n ifi in the range ~ 0. 15 — 0.35. 



V(ho) is a uniform distribution of initial mass-to-flux ratios ho — 
0.1 — 1.0, V(fh c ) is calculated numerically. (The lower limit on ho 
may be as large as ~ 0.6 without changing the high-mass end of the 
predicted CMF; see below. Also, other ho PDFs, e.g., a Gaussian or 
even an inverted Gaussian, give essentially the same CMF, provided 
that they are broad enough to include the range 0.5 — 1.0.) 

In Figure[3] we plot the number of cores N oc V(rh c ) (nor- 
malized to unity) versus m c (normalized to the initial thermal crit- 
ical mass mx,o) for P ex t = 0.1 and r n i,o = 0.162 (dashed 
line), 0.230 (solid line), and 0.325 (dotted line). We also give the 
corresponding a = 1 — d In N/d In m c obtained from a least- 
squares fit to the high-mass end of each curve. For comparison 
with observations, we superimpose data points from SCUBA ob- 
servations of starless cores in O rion at 450 H m an d 850 /im 
dNutter & Ward-Thompson] l2007t) . The width of the mass bins is 
as in Nutter & Ward-Thompson and is shown as a horizontal bar 
for each data point. The vertical error bars denote the V^V count- 
ing uncertainty due to the number of cores in each mass bin. In 
order to match the maximum of the predicted CMF and that of the 
observed CMF, we have taken m,T,o — 2 Mq. 

The excellent agreement between theory and observations is 
evident. A striking feature of the predicted CMF is its lognormal- 
like shape, despite the input uniform distribution of mass-to-flux 
ratios; i.e., one does not need a lognormal density distribution in or- 
der to generate a lognormal-li ke CMF. This is an importa nt distinc- 
tion from past work (e.g., bv lPadoan & Nordlundll20o"2h . Another 
important difference is that in our predicted CMF gravity plays a 
central role; in the Padoan & Nordlund CMF, gravity is irrelevant 
by assumption. 

The predicted CMF has several noteworthy features: 

(i) There is a definite mass turnover. The peak dimensionless 
mass corresponds to cores that had initial mass-to-flux ratios Ho — 
0.6. This result is insensitive to the free parameters P ext and r n i,o- 

(ii) The functional form of the predicted CMF is independent 
of the minimum initial mass-to-flux ratio /io.mm, provided that 
Mo.min < 0.6. The low-mass end of the predicted CMF has a slope 
exactly equal to unity (i.e., a = for m c < 0.3 mT,o). The pre- 
dicted slope of the high-mass end requires only a near-uniform PDF 
of initial mass-to-flux ratios within a factor of ~ 2 smaller than crit- 
ica l, a requirement entirely within observational constraints (e.g., 
see lHeiles & Crutchedl 20051) . 



CMF due to Ambipolar Diffusion L5 



(iii) There is a correlation between the high-mass slope a and 
the maximum core mass m Cjmax /mT,o, which is shown in Fig- 
ure [4] This relation, however, is only indirect: a smaller (greater) 
T n i,o goes hand-in-hand with both a steeper (shallower) high-mass 
slope a (Fig. [3j and a greater maximum core mass m c , m ax/w.T,o 
(Fig.|2j. In other words, molecular cloud cores that form in regions 
with relatively large (small) ionization-equilibrium parameter K, 
(oc ni/nj/ 2 ) should have relatively steeper (shallower) CMFs and 
relatively greater (smaller) maximum dimensionless masses. 



4 SUMMARY AND DISCUSSION 

We have formulated a novel way of predicting the initial mass func- 
tion of molecular-cloud cores (CMF), in which magnetic fields and 
ambipolar diffusion play a central role. The results are in excel- 
lent agreement with observations, provided that the values of the 
two important free parameters remain within their observationally- 
suggested limits. These two parameters, which characterize the par- 
ent cloud at the onset of fragmentation, are: the mass-to-flux ratio 
/xo in units of its critical value for collapse, and the neutral-ion col- 
lision time r n i,o in units of 0.385rfl,o, where Tff,o is the (spherical) 
free-fall time. The results are insensitive to the values of a third free 
parameter (P ox t) for typical molecular clouds, which are gravita- 
tionally bound rather than pressure-confined objects. We find that: 

(i) For typical molecular cloud conditions, the slope a in the 
power-law relation dN/dm c oc is ~ 2.1 — 2.5, with 2.35 
being typical, for core masses m c > 1.5 tot.o- 

(ii) Molecular cloud cores that form in regions with relatively 
large (small) ionization-equilibrium parameter K, (oc rii/nj/ 2 ) 
have relatively steeper (shallower) CMFs and relatively greater 
(smaller) maximum dimensionless masses (in units of the thermal 
critical mass in the parent cloud). 

(iii) There is a mass turnover at m c ~ 0.7 ?71t,o, while for core 
masses m c < 0.3 m-r.o, a = 0. The mass turnover and the low- 
mass slope of the CMF are independent of the free parameter T n i,o- 

For typical molecular cloud conditions (n n .o — 10 4 cm -3 , To ~ 
10 K, P cxt ~ 0.1), the thermal critical mass is mT,o — 4.7 M©. 
Therefore, we expect a = below ~ 1.4 Mq, a mass turnover at 
~ 3 Mq, a ~ 2.1 - 2.5 (with 2.35 being typical) above ~ 7 M Q , 
and an upper cutoff ~ 24 — 47 Mq . 

Some cautions are in order: First, observationally derived core 
masses will be smaller than or equal to core masses derived here, 
depending primarily upon the tracer used (i.e., observations detect 
only a fraction of the total gravitationally-bound core mass). This 
does not affect the overall shape of the CMF, but shifts it toward 
smaller masses. Second, the predicted CMF does not take into con- 
sideration the fact that fragments with greater /iq will evolve toward 
the critical state faster than fragments having smaller fiQ. This im- 
plies that the predicted mass turnover is robust, but observed cores 
with masses below the turnover will most likely be fewer than the 
predicted numbers. Third, the predicted CMF refers to the initial 
core mass function, i.e., it ignores any subsequent fragmentation 
of cores into binary or multiple systems. However, recent observa- 
tions indicating that the functional forms of the stellar IMF and the 
CMF are very similar suggest that this effect may not be important. 

Altogether, then, it presently appears that cosmic magnetism 
and ambipolar diffusion play a crucial role in determining the ini- 
tial core mass function. 
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